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A generalized Bogoliubov model of the Bose gas in the ground state is proposed which properly 
takes into account both the long-range and short-range spatial boson correlations. It concerns 
equilibrium characteristics and operates with in-medium Schrodinger equations for the pair wave 
functions of bosons being the eigenfunctions of the second-order reduced density matrix. The 
approach developed provides reasonable results for a dilute Bose gas with arbitrary strong interaction 
between particles (the 'strong-coupling' case) and comes to the canonical Bogoliubov model in the 
weak-coupling regime. 
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I. INTRODUCTION 



It is well-known that the Bogoliubov model (BM) of 
the weakly interacting Bose gas |l]] is a fundamental of 
the theory of the many-boson systems. The long-range 
spatial correlations of bosons are properly taken into ac- 
count within this model. As to the short-range ones, 
there are situations when the latter need more accurate 
treatment. In particular, one can mention the troubles 
appearing within BM when arbitrary strong repulsion be- 
tween bosons is expressed in a nonintegrable (singular) 
interparticle potential $(r) behaving at small separations 
as l/r m (m > 3) . These troubles are commonly over- 
come by means of using BM with an effective, 'dressed' 
interparticle potential (instead of the 'bare' one, $(?') ) 
that contains all the necessary 'information' concerning 
the short-range boson correlations [^|-[|. 

Although there exist sufficient amount of comprehen- 
sive ways of constructing the effective interaction poten- 
tials ( the pseudopotential method Q , various procedures 
based on summation of the ladder diagrams [[|D ) 1 , it 
looks interesting and promising to realize an alternative 
variant of taking account of the short-range boson cor- 
relations. We mean a generalization of BM which op- 
erates directly with the 'bare' potential $(r) and pro- 
vides a reasonable treatment of the short-range particle 
correlations side by side with the long-range ones. A 
generalization of BM like this is proposed in the present 
Letter. Note that we limit ourselves to the case of the 
zero temperature and consider only equilibrium charac- 
teristics such as the pair correlation function and boson 



momentum distribution. 

The key point of generalizing BM is based on reject- 
ing the usual way of dealing with the Bogoliubov model. 
The investigation presented concerns the second-order 
reduced density matrix (2-matrix). In particular, we op- 
erate with in-medium Schrodinger equations whose solu- 
tions are the eigenfunctions of the 2-matrix, or the pair 
wave functions. As an in-medium interparticle poten- 
tial depends on these functions, so the cited Schrodinger 
equations are nonlinear ones. However, they can be lin- 
earized in the weak-coupling regime as well as for a dilute 
Bose gas even with strong repulsive interaction between 
bosons. The former case corresponds to the canonical 
Bogoliubov model. The latter variant is related to, say, 
the 'strong-coupling' generalization of BM. 

The present Letter is organized as follows. In the sec- 
ond part BM is reconsidered in the framework of the 
2-matrix. The third section concerns the pair wave func- 
tions in the generalized BM. At last, to show a reason- 
able character of the approach proposed, the zero-density 
limit for the Bose gas with strong repulsion between 
bosons is discussed in the fourth part of the paper. 



II. THE BOGOLIUBOV MODEL IN THE LIGHT 
OF THE 2-MATRIX 

Let us consider a homogeneous cold many-body sys- 
tem of N spinless bosons with the volume V and inter- 
particle potential <E>(r) . Note that absence of the spin 
degrees of freedom simplifies the further reasoning with- 
out a loss of generality. 



1 A11 the procedures can be reduced to the ordinary two- 
particle Schrodinger equation with the 'bare' potential. 
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To start our investigation, let us recall the necessary 
definitions. The 2-matrix for the system of interest is 
written as follows: 



f>2 



(ri,r 2 ;ri,r 2 ) = ^ w v ip v (r[, r 2 ) ?/>*(ri,r 2 ) , (1) 



where ^„(ri, r 2 ) are usually called || the pair wave func- 
tions and, physically, w v — 1 (w„ > 0) . The pair 

V 

wave functions for bosons are symmetric with respect to 
the permutation of particles and obey the standard nor- 
malization condition 

/ / |V,(ri,r 2 )| 2 d 3 nd 3 r 2 = l. 
Jv Jv 

The 2-matrix is connected with the pair correlation func- 
tion 

F 2 (r 1 ,r 2 ;r' 1 ,r' 2 ) = ^ + (ri)^ + (r 2 )V(r 2 )V(r' 1 )) (2) 
by the expression Q 

F 2 (t 1 ,t 2 ;v' 1 ,v' 2 ) 



/9 2 (r 1) r 2 ;r 1 ,r 2 ) 



N(N — 1) 



(3) 



Here ip{ r i) denotes the boson field operator. Knowing 
the pair correlation function, one is able to calculate all 
the important thermodynamic quantities . 

The most general structure of the 2-matrix of the equi- 
librium many-body system of spinless bosons is given by 
the following expression ||||: 

. _ pj 



x exp{iq(R' - R)} 



E -yr^r^p^r') exp{zq(R' - R)}, 



(4) 



p q 



where r = ri — r 2 , R = (ri + r 2 )/2 . The quantity 
<Puj,q(r) ■ exp(iqH)/W denotes the wave function of the 
oj— th bound state of the pair of bosons with the total mo- 
mentum hq . Respectively, ip PyCl (r) • exp(«qR)/V stands 
for the wave function of a dissociated state of the pair 
of bosons with the total momentum hq and the momen- 
tum of relative motion hp . For the characteristics N u>q 
and iVp jq we have: N u ^ q is the duplicated number of the 
bound pairs of the w— th species with the total momen- 
tum hq ; -/Vp q is the duplicated number of the dissoci- 
ated pairs with the total momentum hq and the momen- 
tum of relative motion hp . The wave functions Lp u ^ q (v) 
and (^p.q(r) obey the normalization conditions 



lim 



lim i / |^ q (r)| 2 d 3 r = l, 



(5) 



and have the symmetry properties 



<^P,q( r ) = ^P,q(- r ) = ^-P,q( r ) , 

which are a consequence of the Bose statistics. Thus, 
the first term in the right-hand side of (|J) represents the 
sector of the bound pairs; the second one corresponds to 
the dissociated states. Remark that generally speaking, 
one can expect a discrete index to appear in addition to 
q and p for the dissociated states in rather complicated 
situations. However, this does not concern our present 
consideration. So, we have restricted ourselves to the 
summation over q and p in the second term of the right- 
hand side of (Q). 

Comprehensive analysis recently fulfilled in paper ||, 
has demonstrated that, in the thermodynamic limit, the 
correlation function (Q) can be rewritten as 



F 2 (ri,r 2 ; r'i,r' 2 ) = nl<p(r)(p(r') + 



d 3 pd 3 q 



n 
(2tt) 3 



,(p-S).( P+ s )+4 ( p+ a) B ( p _s; 



(^(r V P (r') exp{iq(R' - R)} + F 2 (n, r 2 ; r\, r' 2 ), 



|W 9 (r)rd-V = 1, 



(G) 



where no denotes the density of the condensed parti- 
cles; n(p) = n(p) stands for the distribution of the non- 
condensed bosons over momenta. Note that the Bose- 
Einstein condensation of particles is accompanied by the 
condensation of the particle pairs and, thus, by the ap- 
pearance of the S— functional terms (the off-diagonal 
long-range order) in the pair distribution over momenta 
h p and h q . The first term in the right-hand side of 
(^) is conditioned by presence of a macroscopic number 
of the pairs with q = p = 0. Since they include only 
the condensed bosons, we can call them the condensate- 
condensate pairs. The second term in (|^) corresponds 
to the condensate-supracondensate couples. Besides a 
condensed particle, they also include a noncondensed bo- 
son. At last, F 2 (ri , r 2 ; rj , r 2 ) is the contribution made by 
the supracondensate-supracondensate dissociated states 
of a pair and, maybe, by its bound states. For the wave 
functions of the condensate-condensate and condensate- 
supracondencate couples we have 

(p(r) = l + ^>(r), (^p(r) = V2cos(pr) +V P (r) (p ^ 0), 

(7) 
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where the boundary conditions 

ip(r) —> Q (r —> oo), V'pC 1 ') — * ( r ~~ * 00 ) 



(8) 



take place. At small particle separations the pair wave 
function y p , q (r) is very close to the usual wave function 
of the two-body problem with the relative momcmtum 
p. Therefore, for a singular interparticle potential, when 
<I>(r) cx l/r m (to > 3) at small r, we have y Piq (r) — > as 
r — > . And, hence, the relations 



^(r = 0) = -l, V P (r = 0) = -V2 



(9) 



are fulfilled. 

In the case of a small depletion of the zero momentum 
state (it is of interest in this Letter) we can neglect the 
third term in expression (^): 



16rio 
(2^)3 



F 2 (r 1 ,r a ;r' 1 ,r , 2 )=n§ ^(r)^(r') 

/" d 3 pn(2p)^ P (r)^ p (r')exp{*2p(R' - R)}. 



(10) 



As it is known, there are two physical situations when the 
Bose condensate fraction is expected to be close to 1. One 
of them is related to the weak-coupling regime when a 
small depletion of the zero momentum state results from 
a weak interaction of bosons. The Bogoliubov model is 
an adequate approach of investigating this case. Another 
situation occurs when we deal with a dilute Bose gas with 
an arbitrary strong interaction between particles (singu- 
lar potential). Here the dilution of the system gives rise 
to the small depletion. In this 'strong-coupling' regime 
the short-range correlations play a significant role, which 
is expressed in relations (j^) . On the contrary, the weak- 
coupling case is specified by the inequalities 



|^(r)|«l, |^ P (r)|«l 



(11) 



In particular, the Bogoliubov model is characterized by 
the choice || 

|^(r)|«l, V P (r)=0. (12) 
Expressions ([!(]) and ( [L2"l) allow one to obtain 



F 2 (ri,r 2 ;ri,r 2 ) 



+ 2n n — n a + 



(2tt) 3 



7(r)=ng(l+V(r)+^*(r) 



i(k) exp(ikr) d 3 k , (13) 



where n = N/V and g{r) is the radial distribution func- 
tion. According to the weak-coupling conditions ( |Tl| ) 
and the approximation adopted in (|To|), it is correct to 
neglect the terms of the order of ip(r)(n— no) and [n— no) 2 
111 @. Besides, we may choose the wave function of the 



pair ground state as a real quantity pi: ip(r) = ip*(r) . 
So, expression (13) can be rewritten as 



g(r) = 1 + 2<tp(r) + , n % [ n(k) exp(ikr) d 3 k. (14) 
(2n)- i n J 

Let us show that (|l4|) does represent the result of BM. 
To be convinced of this, we need the equality 

i/)(k) = [ ip(r) exp(-ikr) d 3 r = — (a k a-k) (15) 
J n 

connecting ip(r) with the boson annihilation opera- 
tors 0. At T=0 (we deal with the zero temperature 
case in the present Letter) BM yields jl||| the following 
relations: 



(ikfl-k) 



i-4 



, n(k) 



1-4 



where 



A k = 



no®(k) ( 



and 



E(k) 



h 2 k 2 

E{k)-- n d>(fc) 

2to 



n $(fc) + — I . 

TO 4to z 



(16) 



(17) 



<J)(fc) = / $(r) exp(— ikr) d 3 r. 
Using (p^), ( p"5| ) and ( |l6| ) one is able to arrive at 



(18) 



9{r) = 1 + j^- I exp( i kr) d 3 k . (19) 

(27r) d n J l-Ak 

This relation is exactly the result of the Bogoliubov 
model (see Ref. |(|). 

Concluding this part of the paper, let us take notice 
of an interesting equation following from ( |T5| ) - Qlq ) and 
being important for the reasoning of the next section. It 
is given by 



h 2 k 2 



ip(k) = $(fc) + 2 $(fc) (n(k) + n a i){k) \ , (20) 



where rioip{k) can be replaced by nif)(k) because we 
agreed to neglect the terms of the order of (n — no) ^(r) . 
From and @ it follows that 

^- V 2 <p(r) = *(r) + n J $(|r - y|) (g(y) - l) d 3 y . 

(21) 

This looks like the Schrodinger equation in the Born ap- 
proximation. 



3 



III. PAIR WAVE FUNCTIONS IN THE 
GENERALIZED BOGOLIUBOV MODEL 

As it has been noted above, this Letter addresses the 
generalization of the Bogoliubov model in such a way 
that the short-range correlations should be taken into 
account properly side by side with the long-range ones, 
a small depletion of the zero momentum state being im- 
plied while generalizing. So, in our further investigation 
it is correct to rely on expression ( |l0| ) taken beyond the 
weak-coupling regime introduced by (0). 



To employ approximation (10) for the 2-matrix, one 



needs to determine the wave functions tf(r) and y p (r) 
beyond the weak coupling. To do this, let us consider 
the in-medium two-particle problem: 



H12 ^(ri,r 2 ) = E v Vv(ri,r 2 ) 



(22) 



The Hamiltonian H\ 2 of two bosons placed into the 
medium of similar bosons can be represented as 



Ui + U 2 



(23) 



where Ui (i — 1, 2) stands for the energy of the interac- 
tion of the i—th particle with the medium. Proceeding 
in the spirit of the Thomas-Fermi approach (for details 
see Ref. jL2j) and, thus, neglecting retarding effects, one 
is able to approximate U{ in the form 

Ui = (N-2) /$(|r i -r 3 |Hr 1 ,r 2 ,r 3 )d 3 r3 (i = 1,2), 



(24) 

where w(ri, r 2 , denotes the density of the probability 
of observing the third particle at the point r 3 under the 
condition that the first and second ones are located at ri 
and r 2 . This quantity is connected with the third and 
second reduced density matrices via the relation 

w(ri,r 2 ,r 3 ) = p% 1 (ri,r 2 ;r 1 ,r 2 )p 3 (ri,r 2 ,r 3 ;ri,r 2 ,r 3 ) . 

(25) 

Using the Kirkwood superposition approximation p3[ 
V 3 p 3 (ri,r 2 ,r 3 ;ri,r 2 ,r 3 ) 
^ffdri-raDffdn-raDsC^-ral) , (26) 
one can obtain 

w(ri,r 2 ,r 3 ) ~ i fir(|ri - r 3 |) g(\r 2 - r 3 |) . (27) 
With @ and @ we arrive at 

Ui = U 2 =n I g(\v - y|) $(|r - y|) g(y) d 3 y . (28) 



Thus, equation ([22|) taken with the specifications fl2£ 
and ( p8| ) separates in the usual variables r = ri r 2 and 
R= (n +r 2 )/2, 



ip v (rx,r 2 ) = ip p (r) 



exp(iqR) 



V 



and yields the following relation for the wave function 
Mr): 



V 2 Mr)+$(r)Mr) 



+ 2<p p (r)n J g(\r-y\)<P(\r~y\)g(y)d 3 y = e p tp p (r). 

(29) 

Let us consider equation (p9|). The quantity s p can be 
found taking the limit r — > oo . Using the asymptotic 
relations 

lim g(r) = 1 , lim $(r) = , 

r — >oo r — *oo 

we readily obtain at r — ► oo 

5 (|r-y|)$(|r- y |) ff (y)d 3 y^ | ff(|r-y|)$(|r-y|) d 3 2/. 
Thus, we arrive at 

h 2 P 2 



in 



+ 2n / g(\v-y\)mv-y\)d 6 y (30) 



rather than e p — fi p /m which appears within the or- 
dinary two-body problem. Inserting ( [jO| ) into (p9|), one 
is able to find 

-V 2 ^ p (r) = --^ + $(r)^ p (r) 

+ 2^n^ p (r) J 5 (|r-y|)<i>(|r- y |)(. g (y)-l)d 3 y, 

(31) 

where the equality 

J 5 (|x-y|)$(|x-y|)d 3 y = | g(\z-y\) $(|z-y|) d 3 y 

is used. In ( |3l| ) is the correcting factor which should 
be introduced to compensate oversimplification of treat- 
ing the exchange effects while deriving (p9|). Indeed, the 
exchange between bosons in the pair is taken into con- 
sideration: (fi p (r) — y? p (— r) = (p_ p (r) . However, the 
arguments resulting in ( |29| ) ignore the exchange between 
the particles of the pair and surrounding bosons whose 
influence is considered on the mean-field level, in the 
spirit of the Thomas-Fermi approach. This correcting 
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factor may be, in general, a function of p and some other 
quantities related to the problem: t; ex — £ e a;(p, ■ ■ •) • F° r 
p = equation ([|l]) has to come to ( pl| ) provided re- 
lations (12) are valid, which allows one to find in the 
weak-coupling regime £ e x(p — 0) = 1/2. In this case 
the approximation £ ex — 1/2 can also be employed for 
p / owing to a small mean-absolute value of a boson 
momentum. Moreover, we expect that in the situation 
of the large condensate fraction, the choice £ ex — 1/2 is 
correct beyond the weak coupling too. The main reason 
for this is that the interaction between a couple of par- 
ticles and the medium is weak in both the cases. For 
the potentials with a repulsive core this is due to a small 
density of surrounding bosons. 

Remark, that in contrast to (|2l]), relation ( |3l|) is re- 
duced to the usual Schrodinger equation of the two-body 
problem in the limit n — > 0. The same occurs for r — ► 
in the situation of arbitrary strong repulsion between 
bosons. Indeed, in this case — > oo at r — > . 

Hence, the second term in the right-hand side of j3l]) 
becomes much less than $(r) at small r . This leads to 
conditions ( pj) fu lfilled at any particle density. It is note- 
worthy that (|3l[) with £ ex = 1/2 coincides with one of the 
basic relations of the approach developed in papers 14 1. 

An important peculiarity of equation (^Tj) is that it can 
be used without any divergency in the integral term in 
the case of a singular interparticle potential because 

5 (r) <j>(r) ^ (r-vO). 

Thus, ([H]) reduced to ( pl|) in the weak-coupling approx- 
imation, well answer our purpose of generalizing BM. 

IV. 'STRONG-COUPLING' CASE 

To deal with the generalization of BM based on ( |To| ) 
and ( |3l| ) with £ ex = 1/2, we need one more equa- 
tion. This is because the number of the functions 
g(r), n(k), <p(r), </? P (r) to be determined is larger than 
the number of the equations at our disposal. Within the 
Bogoliubov model for a weakly interacting Bose gas, the 
relation additional to (|lj) and ( pl| ) at zero temperature 
is of the form 



Therefore, in the limit n->0we arrive at 



n 2 if> 2 (k) = n(k) (n(k) + 1) 
(see relations 



(32) 

L5|) - ©). Remark that © follows from 
the canonical character of the well-known Bogoliubov 
transformation . The question now arises if one may 
employ ( p2[) beyond the weak coupling or not. It turns 
out that (|32|) yields quite reasonable results even in the 
case of a dilute Bose gas with strong repulsive interac- 
tion between bosons. To be convinced of this, let us 
consider equations (|l0|), (|3l| ) and ([32]) in the 'strong- 
coupling' regime. From (B2|j it follows that 



n(k) = \ 



l + 4n 2 ^ 2 (k) - 1 



= ntp Q (k) , 



where i}>o(r) obeys equation ( plj ) taken at n 
p = 0: 



(33) 



and 



! (l + ^oW) - (l + Vo(r)) $(r) . (34) 



The relation (|33|) suggests that all the bosons are con- 
densed in the zero-density limit. So, the use of j32] ) 
does not contradict the common expectation concern- 
ing a large condensate fraction in a dilute Bose gas with 
strong repulsive interaction. According to equation (|3~l|), 
for sufficiently low values of p we have 



^p(r) ~ \/2 tp(r) cos(pr) . 



(35) 



Since ( |33| ) is valid at small boson densities n, one can take 
approximation ([55]) to investigate the thermodynamics of 
a dilute Bose gas. Inserting ( pilf ) into ([!(]) we obtain the 
following ansatz: 



l(r) = <p 2 (r)(l + 



(27r) 3 n 



n(k)exp(ikr)d 3 k) , (36) 



where ip(r) is given by equation (|3l]) at p = and 
^ = 1/2: 



h 2 

— V 2 </j(r) = $(r)<p(r) 
m 



:<p(r) J 9(\r-y\)H\r-y\){g(y)-l)d 3 y. (37) 



This ansatz can be used for arbitrary strong repul- 
sion between bosons without any divergency because 
ip(r) $(?') — > at r — > while $(r) — > oo . It is worth 
noting that ansatz (Q) is also good for a weakly inter- 
acting Bose gas. Really, ( |3^ ) is reduced to (14) with the 
assumption IV'MI ^ 1 an( l neglect of the term contain- 
ing the product ijj(r)n(k) . At last, using ( p6| ) and ( |37| ) 
and taking the zero density limit, one can derive 



lim g(r) = ip (r) , 

n— »0 



(38) 



where <po(r) — l + ipo(f) ■ Equality ( pq ) is the well-known 
result for the pair distribution function of the Bose gas 
of strongly interacting particles pi . 
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V. CONCLUSION 



The generalization of the Bogoliubov model of the cold 
Bose gas has been proposed which is based on equations 
(|32|), ( |36| ) and (|37|). They come from the more com- 
plicated set of equations ( |l0|) , ( |3l| ) with t; ex — 1/2 and 
( |32| ) provided the ansatz (|35|) is used. The generaliza- 
tion properly takes into account the short-range boson 
correlations side by side with the long-range ones. The 
proposed approach yields reasonable results in the weak- 
coupling regime as well as in the 'strong-coupling' case 
of a dilute Bose gas with arbitrary intense repulsion. The 
detailed analysis of the latter variant will be fulfilled in 
the forthcoming paper. As it has been noted in the In- 
troduction, the in-medium Schrodinger equations can be 
linearized not only in the weak-coupling approximation. 
This can also be done in the 'strong-coupling' regime. 
In the former case ([ST]) is reduced to equation ( pif ) linear 
in ip(r) . While in the latter situation (|37| ) comes to an 
equation linear in £(r) = ip(r) — (fo (r) due to the obvious 
inequality |C(r)| < \(p (r)\(n -> 0) . 
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